Saddles and softness in simple model liquids 
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We report a numerical study of saddles properties of the potential energy landscape for soft 
spheres with different softness, i.e. different power n of the interparticle repulsive potential. We 
find that saddle-based quantities rescale into master curves once energies and temperatures are scaled 
by mode-coupling temperature Tmct, confirming and generalizing previous findings obtained for 
Lennard- Jones like models. 



The analysis of the potential energy landscape (PEL) 
of model liquids has allowed to clarify many interest- 
ing phenomena of the supercooled liquid regime and the 
slowing down of the dynamics 0, ll H ffll II llli- 
More recently a promising PEL description was obtained 
studying the properties of saddle points J., .10, 11, l3 |. 
In a previous work some of us reported a numeri- 
cal investigation of the PEL for different Lennard-Jones 
(LJ) like model liquids, focusing on the properties of sad- 
dle points visited by the systems at different tempera- 
tures. The main findings of that work were: i) the exis- 
tence of master curves for saddle-based quantities (sad- 
dle order vs. T, energy elevation of saddles from 
underlying minima vs. Us) when temperatures and en- 
ergies are scaled by mode-coupling temperature Tmct 
(we set ks—l)', ii) a nearly constant ratio between ele- 
mentary saddle energy-barriers AE (energy elevation of 
saddles of order 1 from underlying minima) and Tmct- 
AE ~ 10 Tmct] Hi) a quantitative relationship be- 
tween AE and the Arrhenius activation energy AEArr 
(obtained from low-T diffusivity) : AEatt — 2 AE. Al- 
though obtained for different models, the reported uni- 
versality was not too surprising due to the similar shape 
of the repulsive part of the pair potential, in particular 
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all having the same 



dependence. For this reason a 



wider class of models must be analyzed in order to show 
the robustness of the reported universality. 

In this note we extend the class of model liquids un- 
der consideration, analyzing soft spheres with different 
power n of the interparticle repulsive potential (different 
softness) We find that systems interacting with r^" 
potential belong to the same universality class of LJ-like 
potentials, thus pointing towards a common organization 
of saddles in the PEL of disordered systems. 

The investigated systems are 80:20 binary mixtures of 
iV=1000 particles enclosed in a cubic box with periodic 
boundary conditions and interacting through the pair po- 
tential 



(1) 



where Q!,/3 £ A,B, aAA = 1-0, (tab = 0.8, aBB = 0.88, 
eAA = 1-0, eAB = 1-5, ess = 0.5 HI Reduced 
units will be used in the following {uaa for length, caa 



FIG. 1: Saddle order 

Tis SI'S & function of T /Tmct for the 
model systems of this work (soft spheres with different soft- 
ness n) together with the data obtained in a previous work 
for LJ-like models [11^ . Dashed line is a power- law fit, while 
full hue is a REM-like fit [l^. 



for energy, {ma^j^/eAAY^^ for time - m is the mass of 
the particles). The analyzed density was p=1.2. The in- 
vestigated values of the parameter n tuning the softness 
of the interaction were n=6, 8, 12, 18. We note that for 
n < 6 crystallization events prevent the study of the su- 
percooled regime. Following previous works 9] , we stud- 
ied the properties of the saddles of the PEL visited by 
the system during its dynamic evolution at a given tem- 
perature (performed through isothermal molecular dy- 
namics simulations with Nose-Hoover thermostat). The 
temperature range investigated is such that, for each n, 
the diffusivity covers about 4 orders of magnitude. Sad- 
dles are located using minimization procedures (LBFGS 
algorithm as implemented by Liu and Nocedal 17]) on 
the pseudo-potential VF=|Vyp {V is the total potential 
energy). Inherent structures are also located (minimiz- 
ing the true potential V) at each temperature. We then 
collected the properties of saddles (energy Cs and order 
Us, defined as the number of negative eigenvalues of the 
Hessian of V) and inherent structures (energy e^g). 
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FIG. 2: Energy elevation (rescaled by Tmct) of saddles 
from underlying minima (eg — ej^) /Tmct against saddle or- 
der Us, for the systems of this work (soft spheres with differ- 
ent softness n) and for LJ-like models [T3 . The straight line 
has slope 30, corresponding to an elementary energy barrier 
l\E /Tmct=1Q- In the inset the same quantity not rescaled 
by Tmct for the soft sphere models. 

In Fig. 1 we show the saddle order Tig QiS Qb function 
of T /Tmct [Tmct is the mode-coupling temperature es- 
timated from the power law divergence of inverse difFu- 
sivity) for different values of softness n. In the figure 
we report also the same quantity for the LJ-like systems 
analyzed in a previous work All the data collapse 

onto a master curve, even if some deviation is present, 
in particular for n=6 and n=18. In a first approxima- 
tion scaling T by Tmct gives rise to a common behavior 
of the analyzed quantity for the different models. It is 
worth noting that a similar master plot has been obtained 



in Ref. TP] for LJ and CS2 systems. 

Looking at the mean energy elevation eg—e^^ of saddles 
from underlying minima as a function of saddle order Us 
(inset in Fig. 2), one obtains nearly straight lines, with 
different slopes for different values of n. This indicates a 
simple organization of the PEL, with saddles equispaced 
in energy above underlying minima, and allows to de- 
fine an elementary energy barrier elevation Ai? from the 
slopes to' of the curves: Ai? — m'/3 (the factor 3 arises 
from the definition of Ug as the number of negative eigen- 
values of the Hessian normalized to the number of degrees 
of freedom 3iV). From the inset in the figure one observes 
that different softness give rise to different Ai? values. In 
the main panel of Fig. 2 the energy elevation, now scaled 
by Tmct, is plotted vs. n^. The data collapse onto a 
master curve, with a mean slope to = m' /Tmct — 30, 
indicating that the different soft sphere models have a 
similar landscape organization with the same elementary 
energy barrier height when expressed in unit of Tmct- 
AE/Tmct — 10. We note, however, a small correlation 
between the slope and the value of n, with lower values of 
n associated to higher values of the slope. Also reported 
in Fig. 2 are the data obtained in Ref. suggesting a 
wider universality of the relationship AE ~ 10 Tmct- 

A last observation arises from the analysis of Arrhe- 
nius energy barriers AEArr, obtained from Arrhenius fits 
exp{—AEArr/T) of the low-T diffusivity data (see Ref. 
[13). Expressed in unit of Tmct one obtains, for dif- 
ferent softness n, values of AEArr in the range 21 — 24. 
Again this finding is in agreement with the previous ob- 
servation for LJ-like models: AEArr — 2 AE. 

In conclusion, the analysis of soft sphere models with 
different softness confirms the results previously obtained 
for LJ-like models, supporting the hypothesis of an uni- 
versal relation controlling the structure of the PEL. 
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